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CLr ABSTRACT 

2 ■ Context. The origin of low-luminosity Type IIP supernovae is unclear: they have been proposed to originate either from massive (~ 25 Mq) or 
' low-mass (~ 9 Mq) stars. 

$3 Aims. We wish to determine parameters of the low-luminosity Type IIP supernova 2003Z, to estimate a mass-loss rate of the presupemova, and 
1 1 . to recover a progenitor mass. 

Methods. We compute the hydrodynamic models of the supernova to describe the light curves and the observed expansion velocities. The wind 
density of the presupemova is estimated using a thin shell model for the interaction with circumstellar matter. 
1 Results. We estimate an ejecta mass of 14.0 +1.2 Mq, an explosion energy of (2.45 ± 0. 18) x 10 50 erg, a presupemova radius of 229 ± 39 Rq, 
f"*^ ' and a radioactive 56 Ni amount of 0.0063 ± 0.0006 Mq. The upper limit of the wind density parameter in the presupemova vicinity is 10 13 
, gem" 1 , and the mass lost at the red/yellow supergiant stage is < 0.6 Mq assuming the constant mass-loss rate. The estimated progenitor mass 
t-H ' is in the range of 14.4 - 17.4 Mq. The presupemova of SN 2003Z was probably a yellow supergiant at the time of the explosion. 

Conclusions. The progenitor mass of SN 2003Z is lower than those of SN 1987A and SN 1999em, normal Type IIP supernovae, but higher 
I . than the lower limit of stars undergoing a core collapse. We propose an observational test based on the circumstellar interaction to discriminate 
' between the massive (~ 25 Mq) and moderate-mass (~ 16 Mq) scenarios. 
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^ , 1. Introduction Recently an interesting subclass of SNe IIP, so called low- 

■ - - > luminosity SNe IIP, was selected observationally (Pastorello et 

Type II plateau supernovae (SNe IIP) with the plateau of ~ 100 al This family is characterized by luminosities, expan- 

days in the light curve are believed to be an outcome of a core sion velocitieSi and rad ioactive 56 Ni masses which are signifi- 

collapse of the 9 - 25 M stars (e.g., Heger et al. [2003. This cantly lower than mose of normal SNe IIR Two yiews on the 

paradigm assuming the Salpeter mass spectrum suggests that Qrigin of i ow _i uminosity SNe nP have been proposed. Turatto 

about 66% of all SNe IIP should be produced by progenitors, et al haye suggested the origin from massive stars> 

i.e., stars on the main sequence, in the 9 - 15 M Q range. At > 25 Mq> which presumably form black holes and eject a low 

present this general picture remains unconfirmed. It could be amount of 5 6 Ni; alternatively, these SNe might originate from 

verified via the determination of ejecta masses from the hydro- low . mass stars? „ 9 Mq> which m expe cted to eject a low 

dynamic modeling for a sufficiently large sample of SNe IIP. amQunt of 5 6 Ni (Chugai & Utrobin |2000l Kitaura et al. [2006 ). 
Unfortunately, only few SNe IIP have the well-observed light 

curves and spectra needed to reliably reconstruct the basic SN The investigation of the origin of these SNe IIP began 

parameters. It is not therefore surprising that up to now ejecta with some confusion. The point is that SN 1997D, the first 

mass has been determined using the detailed hydrodynamic low-luminosity SN IIP, was detected long after the explosion 

simulations for only SN 1987 A and SN 1999em. It is note- and, therefore, was erroneously cl aimed to possess a short 

worthy that, from the point of view of explosion mechanism, (~ 50 days) plateau (Turatto et al. I19981 I. This, in turn, pro- 

SN 1987A is a normal SN IIP; it has the ejecta mass, the ex- voked a conclusion that SN 1997D originated from a low-mass 

plosion energy, and the amount of ejected 56 Ni comparable to (~ 9 Mq) main-sequence star (Chugai & Utrobin [20001 . The 

those of SN 1999em. subsequent discovery of several low-luminosity SNe IIP with 

a long plateau of ~ 100 days (Pastorello et al. 2004) falsified 

Send offprint requests to: V. Utrobin, e-mail: utrobin@itep . ru this conclusion. Until now there were no attempts to model hy- 
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drodynamically low-luminosity SNe IIP on the basis of new 
observational data. The estimate of the ejecta masses of low- 
luminosity SN 1999br and SN 2003Z was made only using a 
semi-analytical model (Zampieri et al. 120031 Zampieri[2005). 
The semi-analytical model, being a sensible tool for the first 
order estimates, cannot, however, substitute the hydrodynamic 
simulations. 

In this paper we study the well-observed low-luminosity 
Type IIP SN 2003Z (Pastorello 120031 Knop et al. [20071 . 
Our approach is based on the hydrodynamic modeling of the 
light curves and expansion kinematics. This paper is organized 
as follows. The observational data are presented in Sect. [2] 
Section [3] describes the modeling procedure. Results, specifi- 
cally, the SN 2003Z parameters and progenitor mass are pre- 
sented in Sect. [4] The implications of the results are discussed 
in Sect. |5] 

A distance to SN 2003Z of 21.68 Mpc is adopted using 
the Hubble constant Ho = 70 kms 1 Mpc 1 and a recession 
velocity of the host galaxy NGC 2742 v cor = 1518 kms , cor- 
rected for the Local Group infall to the Virgo cluster and taken 
from the Lyon Extragalactic Data base. There are no observa- 
tional signatures of the interstellar absorption in the host galaxy 
(Pastorello 2003) so a total extinction is taken to be equal to the 
Galactic value A B = 0.167 (Schlegel et al. [19981 . 

2. Observational data 

SN 2003Z was discovered very young, and a detection limit 
obtained nine days before the discovery allows us to constrain 
the explosion epoch very well, with an uncertainty of only a 
few days (Boles et al. 2003). We adopt an explosion date of 
SN 2003Z to be JD 2452665. The spectroscopic and photo- 
metric monitoring of SN 2003Z started relatively late, about 3 
weeks after the discovery, and most of data were collected at 
the 3.5 m Telescopio Nazionale Galileo in La Palma (Spain) 
and the 1.82 m Copernico Telescope in Asiago (Italy). The 
plateau phase is quite well sampled (also thanks to the am- 
ateur's data), and a few additional data were collected dur- 
ing the post-plateau phase. All data have been reduced using 
IRAF tasksj. In particular, photometric measurements were 
performed using the point spread function fitting technique. 
The SN magnitudes were then computed with reference to a 
sequence of stars in the SN vicinity, and the photometric zero- 
points of the different nights were obtained from observations 
of several standard fields (Landolt [l9921 l. The detailed descrip- 
tion of the data reduction techniques can be found in Pastorello 
et al. (I2"0"07i All observational data of SN 2003Z, first pre- 
sented in Pastorello (120031 1. will be extensively analyzed in a 
forthcoming paper. 

3. Model overview 

A numerical modeling of SNe IIP exploits radiation hydrody- 
namics in the one-group approximation (Utrobin 2004). When 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



applied to SN 1999em (Utrobin |2007l , the code results in the 
basic SN parameters similar to those recovered in the hydrody- 
namics with the multi-group approach (Baklanov et al.[2P05). 
This suggests that the results of our modeling are not hampered 
by one-group approximation of the radiation transfer. 

The presupernova (pre-SN) structure is set as a non- 
evolutionary red supergiant (RSG) star in hydrostatic equilib- 
rium. The adopted composition of the hydrogen envelope is 
solar. We admit a mixing between the helium core and the hy- 
drogen envelope similar to that in SN 1987A and SN 1999em 
(Utrobin 2005] 120071 1. The model explosion is initiated by a 
supersonic piston applied to the bottom of the stellar envelope 
outside the collapsing 1.4 Mq core. The time scale for the pis- 
ton action is ~ 10~ 3 sec. The explosion energy E is defined as 
a difference between the total energy input and the modulus of 
the binding energy of the envelope outside the collapsing core. 
Since in SNe IIP the radiated energy is small compared to the 
explosion energy, the kinetic energy of the freely expanding 
ejecta is practically equal to the explosion energy. 

The strategy of the search for the optimal model is based 
on the specific dependence of the SN light curve and expan- 
sion velocities on the model parameters (cf. Utrobin 2007 ). The 
56 Ni mass is determined by the bolometric luminosity of the 
radioactive tail, while the plateau duration, the plateau lumi- 
nosity, and the velocity at the photosphere depend on the com- 
bined effects of the ejecta mass, the explosion energy, and the 
pre-SN radius. Specifically, the larger ejecta mass results in the 
longer plateau and the lower plateau luminosity; the larger ex- 
plosion energy leads to the higher plateau luminosity and the 
shorter plateau; the larger pre-SN radius results in the higher 
plateau luminosity and the longer plateau. The extent of the 
56 Ni and helium mixing also affects the light curve, particu- 
larly the shape of the plateau at the end of the photospheric 
epoch (Utrobin 2007). We emphasize the role of the maximal 
velocity of ejecta in the constraining of the model, because the 
density in the outermost layers is sensitive to SN parameters 
(Utrobin & Chugai 2005). Furthermore, the effect of a deceler- 
ation of outer layers due to the ejecta/wind interaction will be 
used to estimate the pre-SN wind density as was done earlier 
for SN 1999em (Chugai et al. l2007T >. 

We rely on the photometric and spectroscopic data reported 
by Pastorello (2003). Although the hydrodynamic model repro- 
duces the V, R, and / light curves satisfactorily, the code with 
one-group radiation transfer is focused on the modeling of the 
bolometric light curve. The observed bolometric light curve of 
SN 2003Z is recovered from BVRI photometry using the black 
body approximation for the SN radiation and the standard fil- 
ter responses. This suggests the determination of the radius of 
the SN photosphere /? p h and the effective temperature T e s ap- 
plying a minimization procedure to the errors between the ob- 
served and calculated filter responses for the adopted distance 
and the reddening. The bolometric luminosity is then calcu- 
lated as L = 4;iJ?p h aTg ff , where cr is the Stephan-Boltzmann 
constant. The data on the photospheric velocity evolution are 
taken from the modeling of the Na I and Ba II lines and from 
Pastorello (2003), while the maximal velocity of ejecta, v max , 
is estimated from the Ha absorption profile. The spectrum on 
day 27.48 (Pastorello |2003) yields v max = 7500 + 200 kms" 1 . 
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Table 1. Physical parameters of hydrodynamic models. 
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The same value we derive from the spectrum on day 9 (Knop 
et al. 2007 ). This indicates a very weak deceleration of the SN 
outer layers between these epochs. 

4. Results 

4.1. Optimal model 

An extensive search in the SN parameter space led us to the op- 
timal model of SN 2003Z (model Opt in Table[TJi with the ejecta 
mass M env = 14 Mq, the explosion energy E - 2.45 x 10 50 
erg, the 56 Ni mass M Ni = 0.0063 Mq, and the pre-SN radius 
Ro = 229 Rq. The initial density distribution in the model 
Opt mimics a heterogeneous structure of the RSG with the 
dense metal/helium core and the extended hydrogen envelope 
(Fig-ID- The helium core is mixed with the hydrogen envelope 



Fig. 2. The mass fraction of hydrogen (solid line), helium (long 
dashed line), heavy elements (short dashed line), and radioac- 
tive 56 Ni (dotted line) in the ejecta of the optimal model. 




l ob, (days) 




Fig. 1. Density distribution with respect to interior mass (a) and 
radius (b) for the optimal pre-SN model. The central core of 1.4 
Mq is omitted. 



Fig. 3. Comparison of the calculated bolometric light curve of 
model Opt (solid line) with the bolometric data of SN 2003Z 
evaluated from the photometric observations of Pastorello 
(2003) (open circles). 



so that the helium abundance drops linearly along the mass co- 
ordinate in the mixing zone (Fig. [2]). In this description the mix- 
ing extent is specified by the mass coordinate of zero hydrogen 
abundance (1.8 Mq). Of note, SN 2003Z similar to previously 
studied SN 1987 A and SN 1999em, requires extended mixing 
between the helium core and hydrogen envelope. 

The bolometric light curve of SN 2003Z is fairly well fitted 
by the optimal model (Fig. [3j. Both empirical and calculated 
V and R light curves (Figs. [4^ and [4J}) show an initial peak 
related to the cooling of hot outer layers of shocked ejecta. The 
amplitude and width of this peak depend on the structure of 
the outermost rarefied layers of the pre-SN envelope. We did 
not try to search for the appropriate density structure of the 
outermost layers to reproduce in detail the initial optical peak 
indicated by the amateur observation^. 

The optimal model reproduces satisfactorily the evolution 
of the photospheric velocity (Fig. |4fc). This evolution is de- 
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Fig. 4. Optimal hydrodynamic model. Panel (a): the calculated 
V light curve (solid lines) compared with the observations of 
SN 2003Z obtained by Pastorello (2003) (open circles) and 
Gary (2003) (open triangles). Panel (b): the same as panel (a) 
but for the R light curves. Panel (c): calculated photospheric 
velocity (solid line) is compared with photospheric velocities 
estimated from modeling the Nal doublet and Ball 6142 A 
profiles (crosses) and absorption minima of the Sc II 6246 A 
line measured by Pastorello (2003) (open circles). Panel (d): 
Ha profile, computed with the time-dependent approach (thick 
solid line), overplotted on the observed profile on day 27.48, as 
obtained by Pastorello (2003 ) (thin solid line). 

termined primarily by the density distribution in the ejecta. 
Remarkably, the power low density distribution in the outer 
layers p oc v~" with n = 9.2 (Fig. [5j is very close to that 
with the power index n = 9 preferred by Knop et al. (|2007) 
on the basis of the analysis of the early spectra. The model 
must also reproduce the maximal velocity v max = 7500 + 200 
km s derived from the blue wing of the Ha absorption in the 
spectrum on day 27.48. This means that the initial velocity v s of 
the boundary shell, which forms because of the radiative damp- 
ing of the shock wave (Grassberg et al. |197U Chevalier[198 1 ), 
should satisfy the inequality v s > v max . In the optimal model 
the boundary dense shell with the mass M s = 1.8 x 10 4 Mq is 
seen in the density distribution at the velocity v s = 7540 km s _1 
(Fig-El which meets the above inequality. It is of note that we 
failed to find a good hydrodynamic model of SN 2003Z with 
the larger boundary velocity, v s > 7540 kms -1 . This consider- 
ation emphasizes the role of the maximal velocity of ejecta in 
the constraining of the hydrodynamic model. 

We calculated the Ha line on day 27.48 in the optimal 
model taking into account time-dependent effects of the hydro- 
gen ionization and excitation (Utrobin & Chugai 2005 ). The 
calculated profile reproduces the blue wing of the absorption 
fairly well (Fig. HJi) although the model is not as good in the 



Fig. 5. The density and the 56 Ni mass fraction as a function of 
the velocity for model Opt at t = 50 days. Dash-dotted line 
is the density distribution fit p oc y~ 9 2 . The arrow marks the 
boundary dense shell at the velocity v s = 7540 km s _1 . 

emission. However, the latter is of minor importance for our 
purpose, the diagnostic of the density in outer layers. Indeed, 
the absorption depth in the blue wing of the Ha line is de- 
termined primarily by the optical depth which, in the time- 
dependent approach, is a function of density in the outer layers. 
The fit of the calculated profile to the observations in the blue 
absorption wing therefore strongly supports the optimal model 
of SN 2003Z. The origin of the disagreement in the emission 
component was studied previously for SN 1987 A (Utrobin & 
Chugai 2005 ) and it was found to be related to our simplified 
description of the ultraviolet spectrum in the range of hv > 3.4 
eV. 

The last column in Table Q] is the total mass of the helium 
core; the value of 3.9 Mq adopted in the optimal model cor- 
responds to a 15 Mq progenitor without rotation (Heger et al. 
2000). The accurate value of the helium core mass is of lit- 
tle significance for the optimal model because the results are 
not sensitive to the helium core in a rather broad mass range. 
A weak sensitivity to the helium core mass has its dark side, 
because it prevents us from recovering the progenitor mass 
on the basis of only the hydrodynamic modeling. In this re- 
gard it would be of interest to consider the case of a pre-SN 
with a massive helium core of 8 Mq (model He8 in Table [TJ 
which is expected for the » 25 Mq progenitor without rotation 
(Heger et al. 2000). The helium core in this model is strongly 
mixed similar to the optimal model (Fig. |6^). Although the 
light curve of model He8 deviates from that of the optimal 
model (Fig. |6j)), the difference is small. We believe that with 
a more sophisticated distribution of the mixed helium and ra- 
dioactive 56 Ni one could reach a closer resemblance to the op- 
timal model. This consideration indicates that a massive pro- 
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Fig. 6. Dependence on the mass of the helium core. Panel (a): 
the chemical composition of model He8 that differs from model 
Opt in the mass fraction of hydrogen (solid line), helium 
(long dashed line), and heavy elements (short dashed line), the 
distribution of radioactive 56 Ni (dotted line) being the same. 
Panel (b): bolometric light curve of model Opt (thick solid line) 
is compared with that of model He8 (dotted line). 

Table 2. Observed parameters of auxiliary hydrodynamic mod- 
els. 
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genitor (25 - 30 Mq) for SN 2003Z cannot be precluded on the 
basis of only the hydrodynamic modeling. Note that the ejecta 
mass for the high-mass scenario should be 14 Mq, so the star 
must lose AM « 10 Mq via the stellar wind, presumably at 
the RSG stage during f lsg « 4 x 10 5 yr (Heger [19981 with a 
high mass-loss rate M = AM/t rsg « 2.5 x 10~ 5 Moyr -1 . This 
provides us a way to verify the high-mass progenitor using the 
estimate of the wind density around the pre-SN in the combi- 
nation with the hydrodynamic model. 

The confidence in the derived model parameters is deter- 
mined by the errors in the distance, the extinction, the plateau 
duration, and the velocity at the photosphere. Unfortunately, 
the errors of the SN 2003Z data are not well defined. Although 
the internal accuracy of the photometric data and the velocity 
measurements is better than 5% (Pastorello 20031 1, the luminos- 
ity error is certainly larger because of uncertainties in the dis- 
tance and the extinction, while the photospheric velocity mea- 
sured from spectral lines and the calculated photospheric veloc- 
ity may refer to physically different layers. We therefore adopt 
rather arbitrarily, the following relative errors: 10% in the bolo- 
metric luminosity, 5% in the photospheric velocity, and 5% in 
the plateau duration. Note that the uncertainty of the 56 Ni mass 
is determined by the error in the bolometric luminosity. At first 
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Fig. 7. Dependence of the bolometric light curve of the optimal 
model (solid line) on the basic parameters: (a) the initial radius, 
model Rms (dotted line) and model Rps (dashed line); (b) the 
ejecta mass, model Mms (dotted line) and model Mps (dashed 
line); (c) the explosion energy, model Ems (dotted line) and 
model Eps (dashed line); (d) the total 56 Ni mass, model Nms 
(dotted line) and model Nps (dashed line). 



glance the error in the luminosity is underestimated, since the 
error in the distance could be as large as 10% (or 20% in lumi- 
nosity) given possible uncertainties in the Hubble constant and 
in the corrected recession velocity of NGC 2742. On the other 
hand, our extensive modeling shows that in order to reproduce 
both the observed maximal ejecta velocity and the light curve, 
one needs to adopt the maximal distance to NGC 2742 (21.68 
Mpc). No much room, therefore, is left for the distance varia- 
tion. This justifies our choice for the luminosity error. 

Using the auxiliary models (Tables Q] and [2] Fig.|7J, we are 
able to translate the adopted errors of the observational param- 
eters into the uncertainties in the pre-SN radius of +39 Rq, the 
ejecta mass of +1.2 Mq, the explosion energy of ±0.18 x 10 
erg, and the 56 Ni mass of ±0.0006 Mq. The ejecta mass thus 
turns out to be in the range 12.8 - 15.2 Mq. Amazingly, us- 
ing the semi-analytical model, Zampieri (2005) found that the 
ejecta of SN 2003Z lies in the range 13-19 Mq which is very 
close to our result. With the adopted errors, the pre-SN mass 
(ejecta plus neutron star) falls into the 14.2 - 16.6 Mq range. 

4.2. Progenitor 

The small difference between the computed boundary shell ve- 
locity of 7540 kms -1 and the observed maximal velocity on 
day 27.48 of 7500 ± 200 km s _1 places a large constraint on the 
pre-SN wind density: the absence of the pronounced ejecta de- 
celeration during the first 27 days implies that the wind should 
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Fig. 8. Evolution of the maximal velocity of the ejecta in the 
optimal model for three values of dimensionless wind density 
parameter a>: 1 (thick line), 0.15 (medium line), and 0.01 (thin 
line). The observed maximal velocity with the corresponding 
uncertainty is shown on day 27.48. 



not be very dense. To find an upper limit of the wind density pa- 
rameter w = M/u (where u is the wind velocity), we compute 
the deceleration of the outer layers of the model ejecta assum- 
ing the thin shell approximation (Chugai et al. 120071 . We intro- 
duce for convenience the dimensionless parameter u> = w/w\ 
where W\ = 6.3 x 10 13 g cm -1 ; the latter corresponding, for ex- 
ample, to the combination of M — 10~ 6 Moyr -1 and u = 10 
km s _1 . The calculated evolution of the maximal velocity of the 
SN 2003Z ejecta for the optimal model is shown in Fig.[8]for 
three cases: ca = 1, 0.15, and 0.01. The case to = 1 predicts 
a too strong deceleration and should be discarded. The wind 
with at = 0.15 is barely consistent with the observed maxi- 
mal velocity within the error. We thus conclude that the density 
parameter for the pre-SN wind of SN 2003Z is a> < 0.15 or 
w < 10 13 gem 1 in CGS units. 

Another look at the wind density issue could be provided 
by the analysis of the contribution of the circumstellar interac- 
tion to the late time luminosity of SN 2003Z. Unfortunately, 
the photometry for this SN is limited by the age of » 200 
days at which only a high upper limit for the density of the 
circumstellar matter can be derived. We can consider, how- 
ever, SN 1997D, another low-luminosity SN IIP, which does 
not show signatures of the circumstellar interaction in the light 
curve till day 400 (Benetti et al. 2001). This suggests the ab- 
sence of a dense circumstellar medium around this SN. We can 
use the circumstellar interaction model (Chugai et al. [2007 ) 
to quantify this statement. Given the similarity of SN 1997D 
and SN 2003Z we use the ejecta mass and kinetic energy of 
SN 2003Z to estimate the upper limit of the wind density pa- 
rameter w ~ 6 x 10 14 gcirT 1 for SN 1997D. It is natural to 
admit that a similar wind is characteristic of the SN 2003Z en- 
vironment as well. This upper limit, however, is of little interest 
in the context of SN 2003Z because it is 60 times higher than 
that obtained above from the early time ejecta deceleration. 

The wind density parameter of SN 2003Z estimated from 
the ejecta deceleration is at least seven times lower than that 



for the normal Type IIP SN 1999em (Chugai et al. [2007l >. 
Assuming that the wind velocities of both pre-SNe are com- 
parable we conclude that the mass-loss rate in the case of 
SN 2003Z is at least seven times lower than that of SN 1999em. 
Given that the mass-loss rate increases with the progenitor 
mass (Neuwenhujzen & de Jager 11990b the difference in the 
wind density implies that the progenitor of SN 2003Z is less 
massive compared to that of SN 1999em. 

In order to convert the pre-SN mass into the progenitor 
mass, requires knowledge of the mass lost via the stellar wind 
on the main sequence and at the RSG stage. The mass lost on 
the main sequence for the ~ 15 Mq star is ~ 0.2 Mq (cf. 
Heger et al. 120001 1. while at the RSG stage the lost mass is 
AM = frsgWM, where f rsg ~ 10 6 yr is the duration of the RSG 
stage for the ~ 15 Mq star (Heger 1998 ). The wind velocity 
is a poorly known parameter. For the RSG stars the velocity is 
in the range of 10 - 40 kms" 1 (Cherchneff & Tielens [T9941 
However, the pre-SN of SN 2003Z with the radius of 229 Rq 
seems to be closer to a yellow supergiant (YSG). Indeed, for the 
15 Mq star luminosity of L * 6 x 10 4 Lq (Heger et al. 120001 . 
the effective temperature of the pre-SN of SN 2003Z would be 
r e jf « 6000 K, i.e., typical for a YSG. The wind velocities esti- 
mated from the YSG observations fall into the range of 30 - 40 
km s" 1 (Lobel et al. 120031 Robberto et al. [19931 . With the max- 
imal velocity u — 40 kms -1 the mass-loss rate immediately 
before the explosion is M < 0.6 x 10 -6 Mq yr" 1 . We assume 
that the mass-loss rate during the whole RSG/YSG stage is con- 
stant on average in agreement with the general wisdom that the 
value of M is mainly determined by the stellar mass and lu- 
minosity (Neuwenhujzen & de Jager 1990). The mass lost at 
the RSG/YSG stage is then < 0.6 Mq. Adopting the mass of 
0.2 Mq to be lost on the main sequence, we conclude that the 
progenitor lost 0.2-0.8 Mq owing to the wind. Combined with 
the pre-SN mass of 15.4 +1.2 Mq this implies the progenitor 
mass of SN 2003Z to be in the range of 14.4 - 17.4 Mq with 
the average value of « 16 Mq. 

5. Discussion and conclusions 

The goal of this study was to recover the basic parameters of 
SN 2003Z and to get an idea about progenitor masses of low- 
luminosity SNe IIP. We estimated the ejecta mass to be 14.0 + 
1 .2 Mq, the explosion energy (2.45 + 0.18) x 10 50 erg, the pre- 
SN radius 229 + 39 Rq, and the 56 Ni mass 0.0063+0.0006 Mq. 
Using the ejecta/wind interaction model, we found the upper 
limit of the density parameter of the pre-SN wind and, assum- 
ing the constant mass-loss rate at the RSG/YSG stage, con- 
strained the progenitor mass by the range of 14.4 - 17.4 Mq. 
The estimate of the wind density thus allows us to avoid the un- 
certainty in the progenitor mass stemming from the weak sen- 
sitivity of the hydrodynamic model to the helium core mass. 

The only normal SN IIP, studied in a similar way to 
SN 2003Z, is SN 1999em. Its pre-SN radius is * 500 Rq, the 
ejecta mass is « 19 Mq, the explosion energy is ss 1.3 x 10 51 
erg, and the 56 Ni mass is a 0.036 Mq (Utrobin 120071 , while 
the progenitor mass is « 22 Mq (Chugai et al. 2007 ). A com- 
parison of SN 2003Z with SN 1999em taken together with the 
fact that the low-luminosity SNe IIP are very similar indicates 
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Table 3. Hydrodynamic models for SN 1987 A, SN 1999em, 
andSN2003Z. 
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that this variety of SNe IIP originates from less massive pro- 
genitors than normal SNe IIP. Parameters of three SNe IIP — 
SN 1987 A (Utrobin|2005]), SN 1999em (Utrobin[2007l Chugai 
et al. 120071 . and SN 2003Z (Table EJ — determined on the 
basis of the similar hydrodynamic modeling, suggest a picture 
in which ordinary SNe IIP originate from massive progenitors 
around 20 Mq, while low-luminosity SNe IIP originate from 
less massive progenitors around ~ 16 Mq, close to the aver- 
age mass of the 9-25 Mq range traditionally associated with 
SNe IIP. If our result for SN 2003Z is correct, then the explo- 
sion energy and the amount of ejected 56 Ni should significantly 
increase for the progenitors between ~ 16 Mq and ~ 20 Mq 
(Fig. [9]). Interestingly, the empirical correlation between the ex- 
plosion energy and the 56 Ni mass for normal SNe IIP has been 
demonstrated by Nadyozhin (2003 ). It should be emphasized 
that these correlations, valid for SNe IIP, may not be applica- 
ble to other SNe II. At least SN 1994W (Type Iln event) was 
found to have a low amount of ejected 56 Ni, ~ 0.015 Mq, but a 
normal explosion energy, * 1.3 x 10 51 erg (Chugai et al. 2004). 

Two alternative conjectures about the origin of low- 
luminosity SNe IIP have been proposed: (1) progenitors are 
massive stars, > 25 Mq (Turatto et al. 11998k (2) these SNe 
originate from low-mass stars, ~ 9 Mq (Chugai & Utrobin 
2000; Kitaura et al. 2006). The present estimate of the progen- 
itor mass of SN 2003Z, 14.4 - 17.4 Mq, is in the disparity with 
both the high-mass and low-mass scenarios. 

Our estimate of the progenitor mass is essentially based on 
the assumption that the mass-loss rate at the RSG/YSG stage 
was constant. In fact, the derived wind density refers to the 
close vicinity of SN 2003Z r < R w = vt ~ 10 16 cm (where 
v ki 7 x 10 8 cm s -1 is the SN expansion velocity in the outer lay- 
ers and t ~ 10 7 sec is the characteristic age of the SN). This lin- 
ear scale with the wind velocity u > 10 kms -1 corresponds to 
the wind history during the latest R w /u < 300 yr before the ex- 
plosion. Most of the RSG/YSG stage, i.e. 4xl0 5 -10 6 yr, is thus 
hidden from our sight. One might suggest that the mass-loss 
rate was substantially more vigorous in the past than immedi- 
ately before the SN explosion. If this questionable possibility 
were the case, the 25 Mq progenitor would lose ~ 10 Mq, and 
the high-mass scenario for SN 2003Z would be saved. Of note, 
the essentially higher wind density of pre-SN in this case sug- 
gests that the dense wind at the large radii could be revealed 
via radio and X-ray observations at the large age, t > 10 yr. 

A remarkable property of all three SNe is strong mix- 
ing between the helium core and the hydrogen envelope in- 
dicated by a low minimal velocity of hydrogen matter (Table 
0). Generally, the model with the unmixed helium core shows 




Fig. 9. Explosion energy (a) and 56 Ni mass (b) versus progeni- 
tor mass for three core-collapse SNe. 

a bump in the light curve at the end of the plateau (Utrobin 
2007); the absence of the bump in the available light curves 
of SNe IIP, in addition to the narrow-topped Ha emission at 
the nebular phase, indicates that substantial mixing between 
the helium core and the hydrogen envelope is a universal phe- 
nomenon for SNe IIP. This mixing could be induced by the con- 
vection during the growth of helium core (Barkat & Wheeler 
119881) or/and during the SN explosion (Kifonidis et al. l2006b . 

It is noteworthy that the radii of all three pre-SNe IIP (Table 
[3]) are notably lower than a radius of typical massive RSG, e.g., 
the radius of * 800 Rq for a Ori (Harper et al. [200TT >. We 
have already shown that the SN 2003Z pre-SN was the YSG 
rather than the RSG. For SN 1999em the 22 Mq pre-SN with 
the luminosity of » 10 5 Lq is characterized by the effective 
temperature of ~ 4700 K, which is intermediate between the 
RSG and YSG values. Interestingly, the pre-SN of SN 2004et, 
a normal SN IIP, was probably a YSG as well (Li et al. l2005l . 
These cases suggest that not only the pre-SN 1987 A, but a pos- 
sibly significant fraction of pre-SNe IIP, experience an excur- 
sion on the Hertzsprung-Russell (HR) diagram from the RSG 
towards the YSG before the explosion. Note, some pre-SNe IIP, 
when becoming the YSG, could get into the instability strip 
and manifest themselves before the explosion as long period 
(P ~ 50 — 100 days) Cepheids. If this is the case, the recov- 
ery of the Cepheid period of a pre-SN on the basis of a set of 
pre-explosion observations could provide us with an additional 
tool for the mass determination from the pulsation period. 

The recovered progenitor mass of SN 2003Z raises a cru- 
cial question: what happens to the 9-15 Mq stars which make 
up a 66% of all the stars from the 9-25 Mq mass range. 
Two conceivable suggestions are: (A) all the 9-15 Mq stars 
explode as low-luminosity, or even fainter, SNe IIP; (B) low- 
luminosity SNe IIP originate only from a narrow range of pro- 
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genitors around ~ 16 Mq, while the rest of the 9-15 Mq stars 
produce different varieties of SNe II. A large sample of SNe II 
which is well observed and studied is certainly needed to dis- 
tinguish between the options A and B. It should be emphasized 
that the initial peak, the end of the plateau, and the radioactive 
tail are of crucial importance for the recovery of reliable SN II 
parameters. 

An alternative approach to determine the pre-SN mass 
exploits archival pre-explosion images of host galaxies and 
stellar evolution tracks on the HR diagram. Interestingly, us- 
ing the archival HST images Van Dyk et al. (2003) and 
Maund & Smartt (120051 ) derived an upper limit of ~ 15 Mq 
and ~ 12 Mq, respectively, for the progenitor of the low- 
luminosity SN 1999br. Given a similarity of known low- 
luminosity SNe IIP these estimates make the high-mass sce- 
nario unlikely. 

It should be emphasized that a test of the mass determina- 
tion method based on the pre-explosion images is needed; for 
instance, the pre-SN light could be partially absorbed in a dusty 
circumstellar envelope. In this regard it would be of top prior- 
ity to determine independently the progenitor mass both from 
the pre-discovery images and from the hydrodynamic simula- 
tions. Among SNe IIP, except for SN 1987A, only SN 2004et 
and SN 2005cs were detected in the pre-discovery images (Li 
et al. 120031 Maund et al. |2005l Li et al. |2006t and have the light 
curves and spectra appropriate for the hydrodynamic modeling 
as well (Sahu et al. 120061 Misra et al. 120071 Pastorello et al. 
2006; Tsvetkov et al. 2006 ). These two challenging cases re- 
quire a detailed study. 
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